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ABSTRACT 

AR Aur is the only eclipsing binary known to contain a HgMn star, making it an 
ideal case for a detailed study of the HgMn phenomenon. HgMn stars are a poorly 
understood class of chemically peculiar stars, which have traditionally been thought 
not to possess significant magnetic fields. However, the recent discovery of line profile 
variability in some HgMn stars, apparently attributable to surface abundance patches, 
has brought this belief into question. In this paper we investigate the chemical abun- 
dances, line profile variability, and magnetic field of the primary and secondary of the 
AR Aur system, using a series of high resolution spectropolarimetric observations. We 
find the primary is indeed a HgMn star, and present the most precise abundances yet 
determined for this star. We find the secondary is a weak Am star, and is possibly 
still on the pre-main sequence. Line profile variability was observed in a range of lines 
in the primary, and is attributed to inhomogeneous surface distributions of some el- 
ements. No magnetic field was detected in any observation of either stars, with an 
upper limit on the longitudinal magnetic field in both stars of 100 G. Modeling of the 
phase-resolve longitudinal field measurements leads to a 3cr upper limit on any dipole 
surface magnetic field of about 400 G. 

Key words: stars: magnetic fields, stars: abundances, stars: chemically peculiar, 
stars: individual: AR Aur 
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1 INTRODUCTION 

HgMn stars are chemically peculiar late B stars, charac- 
terized by strong overabundances of Mn, by up to a thou- 
sand times solar, and strong overabundances of Hg, by up 
to a hundred thousand times solar. Strong overabundances 
of Xe, iron-peak elements, Ga, and underabundances of He 
are also often seen in these stars. These chemical anoma- 
lies are beUeved to be produced by the selective radiative 
levitation and gravitational s ettling of chemical elements in 
the outer stellar layers (e.g. iMichaud. Reeves fc Charlandl 
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1 1974 ). However, despite a significant number of observa- 
tional and theoretical studies, the physical processes respon- 
sible for these strong chemical peculiarities remain poorly 
understood. HgMn stars have been generally thought not 
to p osses strong (dyna mically important) magnetic fields 
(e.g. IShorlin et al.ll2002l ). setting them apart from magnetic 
chemically peculiar Ap and Bp stars. 

Recently, line profile variability has been detected in 
some spectral lines of a number of HgMn stars. This has 
been interpreted a s an inhomogeneous surface distribution of 
specific elements IIRyabchikova, Malanushenko fc AdelmanI 



1 19991 : lAdelman et al.ll2002l : iKochukhov et al.ll2005l ). While 
magnetic chemically peculiar Ap and Bp stars display a wide 
range of surface abundance inhomogeneities, HgMn stars are 
the only type of (apparently) non-magnetic A or B stars in 
which such inhomogeneities have been found. 

In Ap and Bp stars, surface abundance inhomogeneities 
are usually attributed to the magnetic field, possibly working 
together with stellar rotation and magnetically-channelled 
mass loss. This has led some authors to suggest that HgMn 
stars may have hitherto undetected magnetic fields, which 
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coul d give rise to the observed abundance inh omoseneit ies 
(e.g. lAdelman et"ai] I2OO2I : iHubrig et al.ll2006l ) A detailed 
examination of the spotted HgMn star a And bv lWade et al.l 
(I2OOQ) found no evidence for such a magnetic field, with an 
upper limit on the dipole field strength of between 50 and 
100 G. Notwithstanding the apparent absence of magnetic 
fields in a And, the more general presence of magnetic fields 
in variable HgMn stars has not been strongly tested. 

Doppler imaging of sur face abundanc e inh omo- 
geneities has been perfo rmed by lAdelman et al.l ([2003) ^-nd 
iKochukhov et all (|2007| ) for t he HgMn star a And. Surpris- 
ingly iKochukhov et al.l ([20071) find secular evolution of the 
Doppler maps on a time scale of years. This is the only 
known case of surface abundance inhomogeneities evolving 
with time in any A or B star. 

AR Aurigae (17 Aur, HR1728, HD 34364) is a unique 
double-lined spectroscopic binary system with nearly iden- 
tical components (B9V-I-B9.5V). Being the only eclipsing 
binary known to host a HgMn star, it allows for a precise, 
model independent determination of fundamental physical 
parameters. The primary of the system has been long known 
as an HgMn star (Wolff & Preston 1978|) and l ine variability 
in the primary has recently been reported bv iHubrig et al.l 
l|2006r) . This makes the system an ideal candidate for a more 
detailed investigation. 

The two main components of the AR Aur system 
have similar masses of 2.48 Mp, and 2.29 Mp, and an or - 
bital period o f 4.135 days (JNordstrom fc JohansenI [l99J) . 
IChochol et al.l ([19881) discovered a third component to the 
AR A ur system, based on thejight-time effect. According to 
lAlbav rak. Ak & Elmasli (2003|). the third star has a mass of 
0.54 Mn, a separation o f 13 A U, and a period of 23.7 years. 
iNordstrom fc JohansenI ([199J) performed a detailed study 
of the photometric light curves, deriving radii, masses, and 
effective temperatures of the two brighter, eclipsing compo- 
nents of AR Aur, as well as considering the age of the sys- 
tem. They suggest that the secondary is likely a pre-main 
sequence star still contracting towards the zero-age main 
sequence line (ZAMS), while the primary is at the ZAMS. 



Table 1. Table of observations. The exposure time includes the 
number of sub-exposures and the length of a sub-exposure. Peak 
S/N is peak signal-to-noise ratio per CCD pixel. 



iRvabchikoval (| 19981 ') performed an abundance anal- 
ysis of both components of AR Aur, based on the 
equivalent width s from photographic spectra reported by 
iKhokhlova et al.l (|l995l ). and found strong overabundances 
of a wide range of elements in the primary, including 
Mn, Fe, Sr, Y, Pt, and Hg, typical of HgMn stars. In- 
conclusive evidence o f the line profile variability i n AR 
Aur A was given bv pTakeda. Takada fc Kitamural (|l979l ) 
and IZverko. Ziznovskv fc Khokhlova (11997 ) b ased o n low- 
quality photographic spectra. IHubrig et al.l (|2006l ) stud- 
ied AR Aur with modern spectroscopic instrumentation, 
demonstrating line variability for several elements and sug- 
gesting inhomogeneous surface distributions for Sr and Y. 

In this work we search for magnetic fields in AR Aur A 
and B, making this system the second HgMn star with very 
precise upper limits on its magnetic field. Additionally, we 
derive improved chemical abundances for the primary and 
the secondary, confirming the primary as a prominent HgMn 
star, and discovering that the secondary is a weak Am star. 



HJD 


Instrument 


Phase 


Exposure 
Time (s) 


Peak 

S/N 


2454070.87384 


ESPaDOnS 


0.590 


4 X 400 


500 


2454074.91547 


ESPaDOnS 


0.568 


4 X 200 


550 


2454076.79173 


ESPaDOnS 


0.021 


4 X 100 


400 


2454084.53282 


NARVAL 


0.894 


4 X 400 


450 


2454089.39584 


NARVAL 


0.070 


4 X 400 


450 


2454090.36577 


NARVAL 


0.304 


4 X 400 


350 


2454091.39868 


NARVAL 


0.554 


4 X 400 


500 



2 OBSERVATIONS 

Observations were obtained with the Echelle Spectropolari- 
metric Device for the Observation of Stars (ESPaDOnS), a 
high resolution spectropolarimeter located at the Canada 
France Hawau Telescope (CFHT), and with NARVAL, 
a nearly identical instrument mounted on the Telescope 
Bernard Lyot at the Observatoire du Pic du Midi, France. 
Both instruments consist of a bench mounted cross- 
dispersed echelle spectrograph, fiber fed from a Cassegrain 
mounted polarimeter unit. They provide near continuous 
wavelength coverage from 3700 to 10500 A at a resolu- 
tion of i? = 65000. Observations were obtained in spec- 
tropolarimetric mode, providing Stokes V spectra as well 
as Stokes I spectra. The da ta were reduced using Libre- 
ESpRIT (jDonati et al.lll997l ). which performs calibrations 
and optimal spectrum extraction, tailored to the ESPaDOnS 
and NARVAL instruments. 

A total of 7 spectra, 3 from ESPaDOnS and 4 from 
NARVAL, were obtained over a period of roughly one 
month. A summary of the observations is presented in Ta- 
bl eU The orbital p hases are calculated with the ephemeris 
of 'Albav rak et"all (|2003l ). HJD = 2452596.4927 + E x 
4^1346656, corresponding to the primary light minimum. 



3 SPECTRAL DISENTANGLING 

The spectra of AR Aur exhibit a complex pattern of two 
very similar absorption line systems, with velocity separa- 
tion changing from zero to ~200 kms~^ on the time scale 
of two days. In addition, one of the components has weak 
intrinsic line profile variability. In this situation a spectral 
disentangling procedure is essential to separate this effect 
from variable line blending due to orbital motion of the bi- 
nary components and to obtain high-quality average spectra 
for abundance analysis. 

We have developed a direct spectral decomposi- 
tion technique similar to the method described by 
iGonzalez fc Levatd (|200y). Our algorithm operates as fol- 
lows. We start with a set of approximate radial velocities for 
each component and initial guesses of their spectra. Then, 
the contribution of component B is removed, all spectra are 
shifted to the rest frame of component A, interpolated on 
the standard wavelength grid and co-added. This yields a 
new approximation of the primary spectrum. The same pro- 
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cedure is used to update the spectrum of the secondary. 
This sequence of operations is repeated up to convergence. 
As the second major step, we use a least-squares minimiza- 
tion to derive radial velocities and adjust continuum nor- 
malization using low-degree polynomials. Calculation of the 
spectra and radial velocities is alternated until the changes 
of these parameters from one iteration to the next are below 
given thresholds. 

This spectral disentangling was applied to overlapping 
20-70 A segments of AR Aur spectra, avoiding Balmer lines 
and regions affected by the telluric absorption. The final 
radial velocities calculated by averaging results obtained 
for different wavelength regions are presented in Table [21 
Fig.[T]illustrates spectral decomposition for the 4485-4525 A 
interval, containing the variable Pt ll 4514.170 A line of 
AR Aur A. 

Disentangling yields separate average high-quality 
{S/N fs 1000) spectra of components A and B, radial veloc- 
ities and the standard deviation curves, characterizing the 
remaining discrepancy between observations and composite 
model spectra. The standard deviation is examined in the 
rest frame of the primary and secondary separately, which 
allows a straightforward and objective identification of in- 
trinsically variable lines by a coherent excess of standard 
deviation corresponding to a certain absorption feature (see 
Fig. [l|. We found that in all cases the intrinsic spectrum 
variability is associated with the lines of AR Aur A. Thus, 
to study the spectrum variability of the primary in detail, 
in the final step of the disentangling procedure we removed 
from the observations the average spectrum of the secondary 
and corrected the orbital radial velocity shifts. 



4 FUNDAMENTAL PARAMETERS 

Effective temperatures for both components of AR Aur were 
determined based on the spectral energy distribution of the 
binary system. The o bserved spectrophotom etry was taken 
from the catalogue of lAdelman et al.l (11989 and synthetic 
fluxes were calculated using th e LLModels model atmo- 
sphere code IjShulvak et al.ll2004l ) . This code produces plane- 
parallel, line blanketed model atmospheres in local thermo- 
dynamic equilibrium (LTE), using an advanced 'line-by-line' 
method for including the effects of me tal lines. The chem- 
ical abundances of iRvabchikovai (|l998l ) were used as input 
for the model atmosphere calculations, to account for chem- 
ical peculiarities. The best fit combination of models gives 
Teff = 10950±150 K for the primary and T^s = 10350±150 K 
for the secondary. The Taft uncertainties were conservatively 
estimated by eye, based mainly on the fit to the Balmer 
continuum. Fig. [2] illustrates the agreement between the ob- 
served and computed spectrophotometry. The error bars for 
the spectrophotometry are smaller then the symbol size in 
this figure. 

The ratio of radii of the AR Aur system was inferred 
from the ratio of Mg ll 4481 A line equivalent widths, as- 
suming identical Mg abundances. This assumption was veri- 
fied during the abundance analysis presented in Sect. [5] The 
Mg abundances for the two components, based on a wide 
range of Mg lines, are within uncertainty of each other, and 
nearly solar in both stars. We find Rb/Ra = 1.033 ± 0.005, 
which is consistent with the value of 1.020±0.015 determined 
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Figure 2. Comparison of the observed (symbols) and computed 
best fit (solid line) composite spectral energy distribution of AR 
Aur for the parameters Tog (A) = 10950 K, Tcfr(B) = 10350 K, and 
Rb/Ra = 1-033. Error bars on the observations are smaller then 
the symbol size. The computed composite spectral energy distri- 
butions, with the temperature of the primary varied by ±150 K, 
are shown as a dashed lines. 



by iNordstrom fc JohanserJ (|l994l l . Our model predicts the 
central surface brightness ratio Jb/Ja = 0.874 in excellent 
agreement with the observed mean value of 0.870 =b 0.006 in 
the V and y bands (|Nordstrom fc Johansenlll994h . 

Our temperatures of AR Aur A and B are in a gree- 
ment with the values of lNordstrom fc JohanserJ ([199J) but, 
probably, are somewhat more precise. While those authors 
included a coarse correction for chemical peculiarities, based 
mostly on magnetic Ap stars, they did not tailor their mod- 
els to the specifics of the (HgMn) AR Aur system and used 
only the photometric temperature indicators. 

We adopt the log q ( A) = 4 .33 and log g(B) = 4.28 
of INordstrom fc JohansenI ([199J), based on their precise 
masses and radii. As noted by these authors, the observa- 
tion that log (?(A) > log g(B) and Ra < Rb suggests that 
the primary is a ZAMS star while the secondary is still con- 
tracting toward the ZAMS. However, this conclusion relies 
on the assumption of identical Mg abundances in the two 
stars, which may not be fulfilled since both components of 
AR Aur are chemically peculiar. The Hertzsprung-Russell 
diagram positions of the two stars are consistent with this 
conclusion, however they do not rule out both stars being 
at the ZAMS. 

The lack of an accurate spectroscopic orbital solution 
based on precise radial velocity measurements of AR Aur 
has prevented inclusion of this system in the list of funda- 
mental eclipsing binary stars (jTorres. Andersen fc Gimenea 
|201(J ). One can expect that new high- quality spectroscopy 
of the AR Aur system presented by iHubrig et al.l (|2006l 'l 
and in our paper might improve on this situation. A com- 
bined set of radial velocities, determined from our spectra 
and taken from Fig. 2 of Hubrig et al., is given in Ta- 
ble [21 The typical uncertainty of these measurements is 
0.5 kms~^. We modelled the radial velocity curves of both 
components adopting the period from lAlbavrak et al.l (|2003l ) 
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Figure 1. Illustration of the spectral disentangling procedure applied to AR Aur. Symbols represent observations and the thick lines 
show best fitting spectra. Resulting disentangled spectra of the components A and B are shown by the thin lines at the bottom of the 
figure. The standard deviation spectra (shifted upwards and scaled by a factor of 3) diagnose intrinsic line profile variability. An excess 
of deviation at A Ri 4514 A demonstrates the presence of weak variability for the Pt II 4514.170 A line belonging to AR Aur A. 
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Figure 3. The orbital variation of radial velocities for AR Aur A 
and B. Symbols show measurements obtained with UVES (cir- 
cles) and ESPaDOnS/NARVAL (diamonds) for the primary (dark 
symbols) and secondary (light symbols). The lines show fitted ra- 
dial velocity curves for the primary (solid line) and secondary 
(dashed line). 



and assuming a circular orbit. We found it necessary to 
correct our radial velocity measurements by 1.32 kms^^ 
to match the zero-point of Hubrig et al. A small phase 
shift of 0.0024 ± 0.0003 with respect to the ephemeris of 
Albayrak et al. is also required to account for the light- 



time effect caused by the third component. The magnitude 
and sign of this shift a re co nsistent with the pred ictions by 
lAlbavrak et"all (|2003l ) and IChochol et al] l|2006l ). The re- 
sulting fits are presented in Fig. O Our semi-amplitudes, 
Ka = 108.54±0.21 kms-i and Kb = 116.99±0.21 kms'S 
yield a mass ratio of Ma/Mb = 1.078 ± 0.003 and masses 
of the individual components Ma = 2.552 ± 0.008 Mq 
and Mb = 2.367 ± 0.008 Mp for the orbital inc lination of 
i = 88.52° ± 0.06° (|Nordstrom fc Johansen|[l994) . 



5 CHEMICAL ABUNDANCES 

A detailed abundance analysis was preformed for both 
AR Aur A and B, using the disentan gled spectra. The 
ZEEMAN2 spec trum synthesis code (JLandstreetl Il988l : 
IWade et al.l 120011 ) was used, which solves the polarized ra- 
diative transfer equations, assuming LTE. Optimizations to 
the code for negligible magnetic fields values were included, 
and a Levenberg-Marquardt x^ minimization routine was 
used to aid in fitting the observed spectrum. 

Atomic data was e xtracted from the Vienna Atomic 
Line Database (VALD) l|Kupka et al.l 11993 ). using an ex- 
tract STELLAR request, with chemical abundances tailore d 
to the two stars, based on the results of lRyabchikoval (|l998l ). 
The model atmospheres calculated for the fundamental pa- 
rameters derived in Sect. |4] using the LLModels code were 
adopted in the abundance analysis. 

Prior to abundance analysis, the average disentan- 
gled spectra of the primary and secondary were corrected 
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Table 2. Radial velocity measurements of AR Aur. Typical un- 
certainties in the radial velocity are 0.3-0.5 kms~ . Ph ases were 
computed using the ephemeris of lAlbavrak et al.l 1 120031) , given in 
Sect. [2] 



HJD 
-2400000 



Instrument Phase 



Va 
(kms^-'-) 



Vb 
(kms^-'-) 



53671.825 


UVES 


0.078 


-24.86 


82.74 


53700.773 


UVES 


0.079 


-24.06 


83.66 


53725.627 


UVES 


0.091 


-30.86 


91.20 


53713.659 


UVES 


0.196 


-75.14 


137.31 


53730.594 


UVES 


0.292 


-76.74 


140.11 


53722.655 


UVES 


0.371 


-51.09 


112.17 


53632.872 


UVES 


0.658 


119.09 


-70.86 


53728.592 


UVES 


0.807 


128.74 


-81.66 


53724.619 


UVES 


0.846 


116.34 


-69.66 


54070.874 


ESPaDOnS 


0.590 


87.93 


-35.01 


54074.915 


ESPaDOnS 


0.568 


74.64 


-20.73 


54076.792 


ESPaDOnS 


0.021 


11.12 


47.13 


54084.533 


NARVAL 


0.894 


95.25 


-42.15 


54089.396 


NARVAL 


0.070 


-19.09 


80.09 


54090.366 


NARVAL 


0.304 


-73.21 


138.50 


54091.399 


NARVAL 


0.554 


66.22 


-11.81 



for the continuum dilution by the other component (see 
iFolsom et al.ll2008i ) using a theoretical, wavelength depen- 
dent Ught ratio predicted by the adopted LLModels atmo- 
spheres. 

The spectra of both stars were fit for chemical abun- 
dances, as well as usini and microturbulence. Fitting was 
performed for 11 independent regions, between 100 and 
200 A long, ranging from 4170 A to 6200 A. The chemical 
abundances used as free parameters in each spectral window 
were chosen by comparing the observation to the atomic 
data by eye, and looking for lines with corresponding fea- 
tures in the observed spectrum. This was then checked using 
synthetic spectra to ensure a good constraint was available. 
Final best fit values are averages of the results over the 11 
windows, with the window-to-window standard deviation of 
the abundance adopted as the experimental uncertainty. The 
final best fit averages and uncertainties are presented in Ta- 
ble [3l In cases where less than four lines are available the 
uncertainty estimate was made by eye, including potential 
normalization errors, blended lines, and the scatter between 
lines. The number of lines of atomic data that contributed 
significantly to the best fit abundances are given in Table [3] 
in brackets. Additional very weak lines were included to pro- 
vide an accurate spectral synthesis, but such lines contribute 
very little to the fit and the derived abundances, and have 
been excluded from these numbers. 



5.1 Abundances in AR Aur A 

The best fit synthetic spectrum for AR Aur A generally re- 
produces the observed spectrum quite well. Sample fits to 
the spectrum of AR Aur A are shown in Fig. (5] The aver- 
age best fit parameters are shown in Table [S] and chemical 
abundances are plotted relative to solar values in Fig. U The 
effect of hyper-fine splitting was examined for Mn lines and 
generally found to be small. In the few cases where hyperfine 



Table 3. Averaged best fit chemical abundances, tisini and mi- 
crotu rbulence ($) for AR Aur A and B as well as solar abundances 
from lAsplund. Grevesse &: Sauvall ll2005h . Values given in brack- 
ets are the number of lines of atomic data on which the abun- 
dances are based. Abundances are given in units of log ( „^ I . 



AR Aur A 



AR Aur B 



Solar 



23.1 ±0.9 kms" 
< 1 kms-l 



22.9 ±0.7 km s" 
< 1 kms~^ 



He 


-2.09 ±0.08 


(3) 


-0.99 ±0.06 


(4) 


-1.11 


C 


-3.86 ±0.18 


(5) 


-3.76 ±0.11 


(5) 


-3.65 


N 






^ -4.2 


(1) 


-4.26 


O 


-3.52 ±0.08 


(4) 


-3.43 ±0.07 


(4) 


-3.38 


Na 






-5.47 ±0.20 


(2) 


-5.87 


Mg 


-4.80 ±0.10 


(8) 


-4.63 ±0.10 


(12) 


-4.51 


Al 


-7.04 ±0.20 


(1) 


-5.37 ±0.10 


(2) 


-5.67 


Si 


-4.32 ±0.20 


(8) 


-4.34 ±0.12 


(8) 


-4.53 


P 


-5.28 ±0.11 


(4) 






-6.68 


S 


-5.09 ±0.10 


(2) 


-4.51 ±0.15 


(4) 


-4.90 


Ca 


-5.85 ±0.20 


(3) 


-6.05 ±0.13 


(6) 


-5.73 


Sc 


^ -9.5 


(1) 


-9.50 ±0.20 


(6) 


-8.87 


Ti 


-6.50 ±0.14 


(45) 


-7.12 ±0.06 


(31) 


-7.14 


V 






-7.71 ±0.20 


(4) 


-8.04 


Cr 


-5.97 ±0.11 


(32) 


-6.03 ±0.17 


(43) 


-6.40 


Mn 


-5.53 ±0.17 


(43) 


-6.39 ±0.11 


(12) 


-6.65 


Fe 


-4.02 ±0.07 


(356) 


-4.38 ±0.08 


(200) 


-4.59 


Ni 






-5.36 ±0.13 


(10) 


-5.81 


Ga 


^ -7.0 


(1) 






-9.16 


Sr 


-7.72 ±0.20 


(1) 


-8.94 ±0.20 


(1) 


-9.12 


Y 


-7.64 ±0.10 


(23) 


-9.32 ±0.40 


(3) 


-9.83 


Zr 


-7.58 ±0.16 


(26) 


-8.93 ±0.15 


(1) 


-9.48 


Xe 


-5.30 ±0.16 


(4) 


^ -6.0 


(2) 


-9.77 


Ba 


-9.48 ±0.18 


(1) 


-8.71 ±0.09 


(3) 


-9.87 


Ce 


-7.81 ±0.12 


(27) 


^ -8.5 


(3) 


-10.34 


Pr 


-8.60 ±0.08 


(6) 


^ -9.2 


(3) 


-11.46 


Nd 


-7.94 ±0.12 


(21) 


-9.18 ±0.20 


(2) 


-10.59 


Sm 


^ -8.0 


(3) 






-11.03 


Eu 


^ -9.5 


(2) 






-11.52 


Pt 


-6.08 ±0.32 


(3) 


^ -7.2 


(1) 


-10.40 


Hg 


-6.29 ±0.33 


(4) 


^ -7.5 


(1) 


-10.91 


Pb 


-6.9 ±0.4 


(2) 






-10.04 



splitting became important, the line was excluded from our 
fit. 

AR Aur A is clearly an HgMn star, with dramatic over- 
abundances of Hg and Mn. The Fe-peak elements are signif- 
icantly enhanced, as are P, Sr, Y and Zr. Several rare earths, 
Xe, Pt and Pb are also strongly enhanced. C and O are ap- 
proximately solar and Ba appears to be near solar as well. 
He is significantly underabundant, as is Al. Only an upper 
limit on the microturbulence of 1 kms^^ could be derived 
for this star, and abundances here have been calculated as- 
suming no microturbulence. A microturbulence of 1 kms~^ 
would decrease the abundances derived by approximately 
0.05 dex. This upper limit on microturbulence is based on 
the window-to-window scatter of best fit values, when mi- 
croturbulence was included as a free parameter in the x^ fit. 
The majority of spectral windows produced a best fit value 
of kms~^, though some windows produced marginally non- 
zero values. Thus the 'best' overall value of kms~^ was 
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Figure 4. Abundances relative t o solar for AR Aur A (circles) and AR Aur B (squares), averaged over all spectral windows modeled. 
Solar abundances are taken from lAsplund et alJ l|2005l V Points marked with an arrow represent upper limits only. 



adopted, and the fitting procedure was repeated using this 
value. 

The ab undance we derive ge nerally agree with those de- 
termined by iRvabchikoval l| 19981 ). The exceptions to this are 
C, Mn, Sr, Y, an d Pt. For these five elements we find lower 
abundances than iRvabchikoval (|l998l ) by between 0.5 and 1 
dex, which is outside of our uncertainty and any likely un- 
certainty in her results. For Sr, Y, and Pt strong line profile 
variability, discussed below (Sect. [S]), could be responsible 
for some of this discrepancy. 



5.2 Abundances in AR Aur B 

For AR Aur B, we again achieve a very good fit to the ob- 
servations, as presented in Fig. (6] The best fit abundances 
are shown in Table [S] and plotted relative to solar in Fig. [l] 

Based on these abundances, AR Aur B appears to be a 
very weak Am star. The Fe peak elements are, on average, 
enhanced by 0.3 dex, while He, C, and O are all nearly solar. 
Ca and Sc are both significantly underabundant, while Ba 
and Nd are both substantially overabundant. Again, an up- 
per limit on the microturbulence of 1 kms~^ was derived, 
and the abundances are based on a kms~^ model. This 
limit was determined in the same fashion as for AR Aur A. 
The underabundant Ca coupled with the overabundant Fe- 
peak elements, Ba, and Nd strongly suggest that this is a 
weak Am star. 

We find a go od agreement between our abundances and 
those derived by IRvabchikoval l| 19981 ) for the majority of el- 
ements. However we find lower abundances for Ti, Mn, Fe, 
and possibly B a, beyond our unce rtainties and the likely 
uncertainties of IRvabchikoval (|l998h . Due to the higher S/N 
and comparatively advanced spectral disentangling used in 
our observations, we consider our results to be more accu- 
rate. 



6 LINE PROFILE VARIABILITY 

We see clear variability in a number of lines in our disentan- 
gled spectrum of AR Aur A. Line variability was diagnosed 
by examining regions with an abnormally large deviation 
between the observations and the model spectra produced 
by the disentangling procedure. Regions in which a large 
deviation was seen throughout most of a line profile were 
considered significant. Regions of deviation outside of a line, 
or in lines with only a small portion deviating significantly, 
were considered to be spurious. In the spectrum of AR Aur 
A such spurious regions are rare, and are almost always at- 
tributable to the edges of spectral orders or cosmic rays/bad 
pixels. Regions with heavy telluric line contamination were 
avoided entirely. This method is superior to that employed 
in previous studies in that it provides a quantitative measure 
of variability, in the form of the standard deviation spectra, 
rather then a simple estimate by eye. 

Clear variability was seen in lines of Cr, Mn, Sr, Y, Ba, 
Pt, and Hg in the primary (see Fig. [T]). These are some of 
the more strongly overabundant elements in the star, with 
the exception of Ba. Lines with observed variability in the 
primary are listed in Table 3) The lines in which we observe 
variability are generally the stronger lines of the elements 
in question. This suggests that all lines of these elements 
are in fact variable, and the cases in which we observe no 
variability simply fall below our detection threshold. 

The rotation and orbital motio n of the AR Aur compo- 
nents is believed to be synchro nized JNordstrom fc JohansenI 
1 19941 : iKhokhlova et al.(ll995l ). Thus, we presume that the 
line profile variations we observe occur with the period of 
4.135 d. Although we cannot derive the line variability pe- 
riod independently, a good agreement between profiles ob- 
tained at the same orbital phase two weeks apart confirm the 
assumption of synchronous rotation. Since particular vari- 
ability pattern is restricted to specific elements (e.g., the 
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Figure 5. Sample fits to the spectrum of AR Aur A. Points represent the observation and the smooth line represents the best fit 
spectrum. Elements which are major contributors to each line have been labeled. 
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Figure 6. Sample fits to the spectrum of AR Aur B. Points represent the observation and the smooth line represents the best fit 
spectrum. Elements which are major contributors to each line have been labeled. 



variability of Hg and Y lines is clear different), we conclude 
that these spectral changes are produced by an inhomoge- 
neous surface distribution of the elements in question. 



In comparison with iHubrig et al.l l|2006l ). we confirm 
their detection of variability in Y, Pt, Hg, and Sr, but we see 
no evidence for variability in O, Na, Mg, Si, C, Ti, Fe, He, 
Nd, and Zr. Additionally, we see evidence for variability in 
Cr, Mn, and Ba which these authors did not note. An exam- 
ple of weak variations seen in the Mn ll lines is presented in 
Fig. [S] A full comparison of specific lines in which variab ility 
is observed was not possible, since iHubrig et al.l (|2006l ) did 



not publish such a list. However, in the lines they did men- 
tion, we see variability in the Pt ll 4061.7, Hg ll 3983.9, Sr ll 
4077.7, and Y ll 4900.1 lines, but we see no clear variabil- 
ity in the He I 5015 . 7 and Nd lll 4927.5 lines. In the cases 
where IHubrig et al.l l|200a ) report variability that we do not 
observe, it may be that the variability simply falls below 
our detection threshold, though this would imply a very low 
amplitude of variability, or that we observed the system at 
unfavourable phases for detecting variations. However it is 
also possible that some of the variability reported is due to 
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Figure 7. Line profile variability in tiie spectrum of AR Aur A. Tlie spectra are plotted for Y II 3982.5 A, Hg II 3983.9 A, Sr II 4215.5 A, 
Pt II 4514.2 A, and Y II 5087.4 A. Disentangled observations for different pfiases (solid lines) are shifted vertically. The dotted line gives 
the average spectrum. The standard deviation profile, shifted upwards and scaled up by a factor of 3, is shown at the bottom of each 
panel. 
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Figure 8. Weak profile variability in the Mn II 4755.7 and 
4764.7 A lines. The format of this figure is similar to Fig. [7] The 
standard deviation profile is scaled up by a factor of 5. 
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spectral features of the secondary that have not been com- 
pletely removed, or contamination by weak telluric lines. 



7 MAGNETIC FIELDS 

The presence (or absence) of magnetic fields in the photo- 
spheres of AR Aur A/B is diagnosed from the presence of 
circular polarisation within spectral lines in the ESPaDOnS 
and NARVAL spectra, produced as a consequence of the 
longitudinal Zeeman effect. 

No circular polarisation signatures were observed in 
any of the individual spectral lines of the primary or sec- 
ondary. We therefore use d Least Squares Deconvolution 
(LSD) (jDonati et al.lll997n to extract mean Stokes 7, V and 
diagnostic A'^ profiles from the 7 individual observed spec- 
tra, and to thereby increase the sensitivity of the magnetic 
diagnosis. 

The LSD analysis was performed using the original re- 
duced spectra representing the combined system. Spectral 
disentangling was not possible for the Stokes V spectra since 
the disentangling technique relies on both detectable lines, 
which we do not have in Stokes V , and lines remaining 



Figure 9. LSD Stokes I (lower) and V (upper) profiles extracted 
from the 7 observations of the AR Aur system. The primary and 
secondary components are distinguished by their slightly different 
line depths. The Stokes V profile amplitudes are scaled up by 25 
times with respect to / for clarity. 



roughly constant, which would not be the case in Stokes V 
if they were detected. The line mask used for the LSD pro- 
cess was constructed using extract stellar requests from 
VALD, for the atmospheric parameters and abundances of 
AR Aur A derived in Sect. U and [S] The mask was filtered to 
include only lines with predicted unbroadened depths great 
than 10% of the continuum, and ultimately contained 1168 
lines. As the atmospheric parameters and Fe-peak element 
abundances of AR Aur A and B are similar, and because 
such lines represent the majority of the lines used in the 
mask, the use of the primary abundances for the LSD ex- 
traction was not expected to be a significant source of error. 
This was confirmed by experiment. The Stokes / and V LSD 
profiles are illustrated in Fig. |5] 

The LSD profiles are characterised by signal-to-noise 
ratios (S/N; per 1.8 km/s pixel) ranging from about 7000 to 
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Table 5. Longitudinal magnetic field measurements of AR Aur. Columns report the UT date of observation, signal-to-noise ratio per 1.8 
km/s LSD pixel in the extracted LSD A'^ profile. Then, for both the primary and secondary respectively; Stokes V detection diagnosis 
(ND=no detection), longitudinal field measured using V and N across the mean line in units of Gauss. On 6 Doc 2006, the profiles of AR 
Aur A and B were blended together. The longitudinal field measured from the combined profiles was 57 ± 33 G (Stokes V) and 18 ± 33 G 
(diagnostic null). 









AR Aur 


A 




AR Aur 


B 


Date 


LSD S/N 


Det 


V Be (G) 


N Bi (G) 


Det 


V Be (G) 


N Bi (G) 


30 Nov 06 


10800 


ND 


-29 ± 26 


-24 ± 26 


ND 


50 ±43 


-77 ± 43 


04 Dec 06 


14200 


ND 


29 ±19 


-31 ± 19 


ND 


-24 ±29 


-10 ±29 


14 Dec 06 


10500 


ND 


-10 ±31 


-3 ±31 


ND 


-31 ± 36 


-19 ±36 


19 Dec 06 


10500 


ND 


2 ±34 


-44 ± 34 


ND 


33 ±48 


27 ±49 


20 Dec 06 


7100 


ND 


9 ±42 


32 ±43 


ND 


12 ±69 


15 ±70 


21 Dec 06 


13400 


ND 


-14 ±23 


-14 ±23 


ND 


-26 ± 36 


-23 ± 36 



about 14000, corresponding to a net S/N gain of about 20 
times. No signal was detected in the V or A'' LSD spectra 
of any of the observa tions, based on the detection criteria 
of iDonati et al.l ij 19971 ). The longitudinal magne tic field was 
derive d from each observation using Eq. (1) of I Wade et al.l 
l|2000f ). To derive the longitudinal field of each of the com- 
ponents of AR Aur, we integrated across the line profiles of 
each star at phases where those profiles were separated in 
velocity and unblended. The integration range was carefully 
established by eye for each profile; the full integration range 
was typically 30 ± 1 km/s on either side of the profile center- 
of-gravity. No statistically-significant (i.e. greater than Sa) 
longitudinal field was detected in either star, with typical la 
uncertainties of 30-40 G. The results of the magnetic anal- 
ysis are reported in Table \^ We point out that the inferred 
value of the longitudinal field of a component of a binary 
system is unaffected by the continuum of the companion. 
This is because it is the continuum normalised Stokes / and 
V parameters {I /Ic and V/Ic, respectively) that appear in 
the denominator and numerator of Wade et al.'s Eq. (1). 
Nevertheless, the inferred uncertainty, derived from photon- 
noise error bars propagated through the LSD procedure, is 
larger than that for the single-star case. 

From these basic magnetic measurements, we conclude 
that neither component hosts a magnetic field with a disc- 
averaged longitudinal component larger than about 100 G. 
If either star hosted a field with a predominantly dipolar 
topology (like those of the magnetic Ap stars), the individual 
measurements indicate that the polar surface strength of 
that field is roughly constrained to be weaker than about 
300 G. 

Because the stellar components of the AR Aur system 
are likely to be tidally locked with their rotation, we can 
phase the longitudinal field measurements and LSD profiles 
according to the orbital (and therefore presumably rota- 
tional) ephemeris give in Sect. 2. We point out that the 
orbital period of 4.1346656 days, interpreted as a rotational 
period, is in excellent agreement with the radii and v sin is of 
the components, assuming that the inclination of the stellar 
rotational axes is equal to the orbital inclination. 

We have phased the longitudinal field measurements, 
and examined their agreement with synthetic longitudinal 
field curves corresponding to a large grid of dipole surface 
magnetic field models. While three of the measurements oc- 



cur at approximately the same phase, the dataset samples 
the entire rotational cycle at approximately quarter-cycle 
intervals. This sparse but relatively uniform sampling is im- 
portant for constraining the field intensity using longitu- 
dinal field measurements. We assume the stellar rotation 
axis inclination i ~ 88.5°, in agreement with the orbital in- 
clination and consistent with our assumption that the sys- 
tem is tidally locked. We find for perpendicular (/3 ~ 90°) 
dipoles, the 6 measurements constrain the polar strength of 
any dipole magnetic field in the primary to below 225 G, at 
3(7 confidence. For intermediate field obliquities (/3 ~ 45°), 
the upper limit on any dipole field present is about 375 G. 
For the secondary, the analogous upper limits are similar: 
about 275 G and 400 G, respectively. 

The lack of detection of any signature in the Stokes V 
profiles also constrains the presence of more complex mag- 
netic fields in the photospheres of AR Aur A/B. However, 
to place quantitative upper limits on the presence of such 
fields requires an ad hoc assumption about their structure 
and intensity. Such an investigation is beyond the scope of 
this paper. 

In summary, our observations provide no evidence for 
the presence of magnetic fields in either of the components 
of AR Aur. If oblique dipolar fields are present in either 
star, the polar strength of those fields are constrained to be 
weaker than a few hundred G at 3cr confidence. As fields 
of inten sity similar to our upper limi t (e.g. 300-400 G in 
e UMa; iBohlender fc Landstreetlll99Gl ) are capable of pro- 
ducing strong chemical nonuniformities in the atmospheres 
of Ap stars, we are not able to fully rule out a magnetic 
origin of these features. Nevert heless, such weak fie lds in Ap 
stars appear to be quite rare - lAuriere et al.l (|2007l ) identify 
only 3 stars out of a modeled sample of 24 stars with fields 
that are inferred to be similarly weak. Therefore, given the 
rarity of such fields, and the high (3a) confidence of our null 
result, we conclude that it is highly probable that the chem- 
ical nonuniformities inferred to exist in the atmosphere of 
AR Aur A are not of magnetic origin, and are produced by 
some other phenomenon. 



8 CONCLUSIONS 

We find AR Aur A to be a st rongly peculiar HgMn 
star, confirming previous results IjKhokhlova et al.l Il995l : 
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Table 4. Variable lines observed in AR Aur A. If a line is blended 
significantly multiple components are listed. Lines marked by a 
question mark show only very weak variability, near the noise 
level. 



Ion Wavelength Ion Wavelength Weak 

(A) (A) 



Yii 


3950.352 










Yii 


3982.594 










Hgii 


3983.93f 










Pt II 


4046.449 


+ 


Hgi 


4046.559 




Pt II 


4061.659 










Sr II 


4077.709 










Yii 


4124.907 








? 


Sr II 


4161.792 










Yii 


4177.529 










Yii 


4204.695 








? 


Mn II 


4206.367 








7 


Sr II 


4215.519 










Yii 


4235.729 










Mn II 


4242.333 


+ 


Cm 


4242.364 


? 


Sr II 


4305.443 










Yii 


4309.631 










Mn II 


4326.639 








7 


Hgi 


4358.323 


+ 


Yii 


4358.728 


? 


Yii 


4374.935 










Yii 


4398.013 










Yii 


4422.591 










Mn II 


4478.637 








? 


Pt II 


4514.170 










Ban 


4554.029 








7 


Cr II 


4554.988 








7 


Cr II 


4558.650 








? 


Cr II 


4618.803 








? 


Cr II 


4634.070 








7 


Yii 


4682.324 










Mn II 


4755.727 








? 


Mn II 


4764.728 








? 


Yii 


4786.580 








7 


Yii 


4823.304 


+ 


Mm 


4823.524 




Cr II 


4876.399 








? 


Yii 


4883.684 










Yii 


4900.120 










Yii 


5087.416 










Yii 


5123.211 


+ 


Mn II 


5123.327 




Yii 


5200.406 










Yii 


5205.724 










Cr II 


5237.329 








7 


Mn II 


5297.028 


+ 


Mn II 


5297.056 


7 


Mn II 


5299.330 


+ 


Mn II 


5299.386 


7 


Yii 


5402.774 










Yii 


5497.408 










Yii 


5509.895 










Yii 


5662.925 











iRvabchikoval 1 1998f ) . AR Aur B shows weak Am peculiar- 
ities, particularly modest overabundances of Fe-peak ele- 
ments and Nd, and modest underabundances of Ca and Sc. 
Thus AR Aur is a binary system with two different chemi- 
cally peculiar stars. This is striking in light of the very simi- 
lar temperatures of the two stars. It is possible that the small 
temperature difference represents a sharp dividing line be- 
tween Am and HgMn stars. Alternately, if AR Aur B is still 



on the pre-main sequence, the difference in chemical abun- 
dances may be a result of the different evolutionary status 
of the stars. 

We confirm the result of iNordstrom fc JohanserJ (119941 ) 
that the cooler, lower mass, secondary has a larger radius 
than the primary, assuming that both components have the 
same Mg abundance. Under this assumption, we agree with 
their conclusion that this implies the secondary is still con- 
tracting towards the ZAMS, while the primary has likely al- 
ready arrived on the ZAMS. Thus both components of AR 
Aur are likely some of the youngest chemically peculiar stars 
known. AR Aur B may be the first pre-main sequence Am 
star discovered. This result, if correct, implies that chemical 
peculiarities must arise quickly in both HgMn stars and in 
Am stars. The atomic diffusion process thought to drive the 
observed chemical peculiarities must be fairly efficient for 
this to occur. 

Despite a careful high S/N search, we find no evidence 
for a magnetic field in AR Aur A or B. The LSD profiles for 
our observations, which are sensitive to simpler field geome- 
tries with a zero net longitudinal component, as well as more 
complex field geometries, show no indication of any signal in 
either star. Our longitudinal field measurements place a Sc- 
upper limit of 100 G on the longitudinal fields of both stars 
over a range of rotational phases. This constrains a dipole 
field to be less than about 400 G at the pole in either star. 
Thus we conclude that neither AR Aur A nor B is likely to 
h ave any signi fi cant magnetic field. This matches the result 
of lWade et ai] (|2006l ') for the HgMn star a And. 

AR Aur represents a second case in which line profile 
variability and surface abundance inhomogeneities exist in 
the absence of any strong magnetic field. This substantially 
strengthens the conclusion that, whatever the mechanism 
producing these inhomogeneities, a magnetic field is not re- 
quired. 
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